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Cavitation and thermal photon production in relativistic heavy ion colhsions 

Jitesh R. Bhattjj Hiranmaya Mishrajj and V. Sreekanttj^ 
Physical Research Laboratory, Ahmedabad 380009, India 

We investigate the thermal photon production-rates using one dimensional boost-invariant second 
order relativistic hydrodynamics to find proper time evolution of the energy density and the temper- 
ature. The effect of bulk- viscosity and non-ideal equation of state are taken into account in a manner 
consistent with recent lattice QCD estimates. It is shown that the non-ideal gas equation of state 
i.e e — 3 P 7^ behaviour of the expanding plasma, which is important near the phase-transition 
point, can significantly slow down the hydrodynamic expansion and thereby increase the photon 
production-rates. Inclusion of the bulk viscosity may also have similar effect on the hydrodynamic 
evolution. However the effect of bulk viscosity is shown to be significantly lower than the non-ideal 
gas equation of state. We also analyze the interesting phenomenon of bulk viscosity induced cav- 
itation making the hydrodynamical description invalid. It is shown that ignoring the cavitiation 
phenomenon can lead to a very significant over estimation of the photon fiux. It is argued that this 
feature could be relevant in studying signature of cavitation in relativistic heavy ion collisions. 



I. INTRODUCTION 

Thermal photons emitted from the hot fireball created 
in relativistic heavy-ion collisions is a promising tool for 
providing a signature of quark-gluon plasma [iHB] (see 
[7H5] for recent reviews). Since they participate only 
in electromagnetic interactions, they have a larger mean 
free path compared to the transverse size of the hot and 
dense matter created in nuclear collisions [10]. There- 
fore these photons were proposed to verify the existence 
of the QGP phase [TTl [T2] . Spectra of thermal photons 
depend upon the fireball temperature and they can be 
calculated from the scattering cross-section of the pro- 
cesses like qq — > (77, hremsstrahlung etc. Time evolu- 
tion of the temperature can be calculated using hydro- 
dynamics with appropriate initial conditions. Thus the 
spectra depend upon the equation of state (EoS) of the 
medium and they may be useful in finding a signature 
of the quark-gluon plasma [T^l - fTF| . Recently the thermal 
photons are proposed as a tool to measure the shear vis- 
cosity of the strongly interacting matter produced in the 
coUisionsdZmH]. 

Understanding shear viscosity of QGP is one of the 
most intriguing aspects of the experiments at Relativistic 
Heavy Ion Collider (RHIC). Analysis of the experimental 
data collected from RHIC show that the strongly coupled 
matter produced in the collisions is not too much above 
the phase transition temperature Tc and it may have ex- 
tremely small value of shear viscosity rj. In fact the ratio 
of the shear viscosity 77 to the entropy density s i.e. 77/5 
is around l/47r which is the smallest for any known liquid 
in the nature[Tn]. In fact the arguments based on AdS- 



CFT suggest that the values of rj/s can not become lower 
than l/47r. This is now known as Kovtun-Son-Starinets 
or 'KSS- bound' [20]. Thus the quark-gluon plasma pro- 
duced in RHIC experiments is believed to be in a form 
of the most perfect liquid [2 IJ. No wonder ideal hydro- 
dynamic appears to be best description of such matter 
as suggested by comparison between the experimental 
data[22, and the calculations done using second-order rel- 
ativistic hydrodynamics 
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However there remains uncertainties in understanding 
the application and validity of the hydrodynamical pro- 
cedure in the relativistic heavy-ion collision experiments. 
It is only very recently realized that the effect of bulk 
viscosity can bring complications in the hydrodynamical 
description of the heavy-ion collisions. Generally it was 
believed that the bulk viscosity does not play a signifi- 
cant role in the hydrodynamics of relativistic heavy-ion 
collisions. It was argued that that since C, scaled like 
e — 3P at very high energy the bulk viscosity may not 
play a significant role because the matter might be follow- 
ing the ideal gas type equation of state[3T]. But during 
its course of expansion the fireball temperature can ap- 
proach values close to Tc. Recent lattice QCD results 
may not have ideal gas EoS and the ratio C/s show a 
strong peak around T^ [Ml 133] ■ The bulk viscosity con- 
tribution in this regime can be much larger than that 
of the the shear viscosity. Recently the role of bulk vis- 
cosity in heating and expansion of the fireball was ana- 
lyzed using one dimensional hydrodynamics 34J . Another 
complication that bulk viscosity brings in hydrodynamics 
of heavy- ion collisions is phenomenon of cavitation[35j. 
Cavitation arises when the fiuid pressure becomes smaller 
than the vapour pressure. Since the bulk viscosity (and 
also shear viscosity) contributes to the pressure gradient 
with a negative contribution, it may be possible for the 
effective fluid pressure to become zero. Once the cavita- 
tion sets in the hydrodynamical description breaks down. 



It was shown in Ref . |35| that cavitation may happen in 
RHIC experiments when the effect of bulk viscosity is in- 
cluded in manner consistent with the lattice results. It 
was shown that the cavitation may significantly reduce 
the time of hydrodynamical evolution. 

One of the main objectives of this paper is to study 
the photon spectra with the effect of bulk viscosity and 
cavitation. Finite ( effect can either significantly reduce 
the time for the hydrodynamical evolution (by onset of 
cavitation) or it can increase the time by which the sys- 
tem reaches TJ Moreover the non-ideal gas EoS can also 
significantly influence the hydrodynamics (see below) . In 
what follows, we use equations of relativistic second order 
hydrodynamics to incorporate the effects of finite viscos- 
ity. We take the value of C/s same as that in Ref. [55] 
and keep 77/s — l/47r. Further we use one dimensional 
boost invariant hydrodynamics in the same spirit as in 
Refs. [331 I3S] . One of the limitations of this approach is 
that the effects of transverse flow cannot be incorporated. 
As the boost-invariant hydrodynamics is known to lead 
to underestimation of the effects of bulk viscosity [53] 5 we 
believe that our study of the photon spectra will provide 
a conservative estimate of the effect. 



II. FORMALISM 

A. Viscous Hydrodynamics 

We represent the energy momentum tensor of the dis- 
sipative QGP formed in high energy nuclear collisions as 



T'"' ^eu'^u'' -PA^"" +11"'' 



(1) 



where e, P and u'^ are the energy density, pressure and 
four velocity of the fluid element respectively. The oper- 
ator A'^'^ = g"" — u" u" acts as a projection perpendic- 
ular to four velocity. The viscous contributions to T^" 
are represented by 



W = TT^" - I\^"' n 



(2) 



where tt^'', the traceless part of 11^'^; gives the contribu- 
tion of shear viscosity and 11 gives the bulk contribution. 
The corresponding equations of motion arc given by, 

De+{e + P)e-WV(^u^)=Q (3) 

(e -I- P) Du'^ - V"P + A„^ d^Ji"" ^ (4) 

where D = u'^df^, 9 = df^u", Vq — A^ad" and 
-^{fi By) = \ [A^ By + Ay Bfj] gives the symmetrization. 

We employ Bjorken's prescriptioniSO to describe the 
one dimensional boost invariant expanding flow, were we 
use the convenient parametrization of the coordinates us- 
ing the proper time r — \/t'^ — z^ and space-time rapidity 



V = \ ^"[f^f ]; t — T coshy and z = r sinhy. Then the 
four velocity is given by. 



u" — (coshy, 0, 0, sinhy). 



(5) 



We note that with this transformation of the coordinates, 
D = ^ and 6* = l/r. 

Form of the energy momentum tensor in the local rest 
frame of the flreball is then given bv |37f[40] 



rp^V 



£000 
Pl 
Pl 
P, 



(6) 



where the effective pressure of the expanding fluid in the 
transverse and longitudinal directions are respectively 
given by 



P± 

Pr 



P 
P 



n + l<i> 



(7) 



Here $ and 11 are the non-equilibrium contributions to 
the equilibrium pressure P coming from shear and bulk 
viscosities. Respecting the symmetries in the transverse 
directions the traceless shear tensor has the form tt*^ = 
diag($/2,$/2,-$). 

In the first order Navier-Stokes dissipative hydrody- 
namics 



n - ~c,d^,u" 



and tt'''' = 277V ^''u''' 
with ^, ry > and ^ i^^Uy) — 2V(^' 



(8) 



|A„^V„u". So 



for first order theories with Bjorken flow we have 

C 4?? 

n = -^ and *=7-^. 

r 6T 



(9) 



The Navier-Stokes hydrodynamics is known to have in- 
stabilities and acausal behaviours [311 132]; second order 
theories removes such unphysical artifacts. 

We use causal dissipative second order hydrodynamics 
of Isreal-Strewart [43] to study the expanding plasma in 
the flreball. In this theory we have evolution equations 
for n and $ governed by their relaxation times rn and 
Ttt. We refer [44l|45] for more details on the recent devel- 
opments in the theory and its application to relativistic 
heavy ion collisions. 

Under these assumptions, the set of equations (i.e., 
equation of motion ([3| and relaxation equations for vis- 
cous terms) dictating the longitudinal expansion of the 
medium are given by[3ni |3H HS] 



de 

dT ~ 


-\{e + P + 


dT 


$ 2 1 

T, + 3 hT 
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n 1 


dT 


Tn (ioT ' 
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(10) 

(11) 

(12) 



where <& = tt''" — tt^^. The terms in the square bracket 



in Equation! 11 1 are needed for the conformaHty of the 



theory ^47j. The coefficients (3o and 132 are related with 
the relaxation time by 



Tn = C /3o , Ttt = 2?7 /32 



(13) 



We use the J^ — A supersymmetric Yang-Mills theory 
expressions for r^ and Ai [47H49] : 



2 - hi 2 
2ttT 



and 



Ai 



27rT 



(14) 



(15) 



We set T7r(r) = rn(T) as we don't have any rehable pre- 
diction for TnjM]. 

In order to close the hydrodynamical evolution equa- 
tions (10 - 12 1 we need to supply the equation of state. 



B. Equation of state, C,j s and 77/s 

We are interested in the effect of bulk viscosity on the 
hydrodynamical evolution of the plasma and recent stud- 
ies show that near the critical temperature T^ effect of 
bulk viscosity becomes important [50i[5T]. We use the re- 
cent lattice QCD result of A. Bazavov et aZ.[31I for equi- 
librium equation of state (EoS) {non-ideal: e — 3P ^ 0). 
Parametrised form of their result for trace anomaly is 
given by 



£-3P 

2^4 



= 1 



1 



exp 



T-ci 

C2 



2^2 



d4 

y4 



(16) 
where values of the coefficients are ^2 = 0.24 GeV^, ^4 = 
0.0054 GeV, ci = 0.2073 GeV, and ca = 0.0172 GeV[i5]. 
Their calculations predict a cross over from QGP to 
hadron gas around .200-. 180 GeV. We take critical tem- 
perature as .190 GeV throughout the analysis. The func- 
tional form of the pressure is given by [32] 



PjT) PjTo) 



n 



dT 



To 



,e ~3P 

2^/5 



(17) 



with To = 50 MeV and P{Tq) = 



From Equations ( 16 ) and ( 17 1 we get e and P in terms 
of T. 

We rely upon the lattice QCD calculation results for 
determining C/s. We use the result of Meyer [3ff, which 
indicate the existence a peak of C/s near Tc, however 
the height and width of this curve are not well under- 
stood. We follow parametrization of Meyer's result from 
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FIG. 1. (e - W)/T'^,C/s (and ri/s = 1/A-k) as functions 
of temperature T. One can see around critical temperature 
{Tc = .190 GeV) C, ^ ri and departure of equation of state 
from ideal case is large. 
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FIG. 2. Various bulk viscosity scenarios by changing the 
width of the curve through the parameter AT. 



Ref . [55] ■ given by 



C (T,-T 

- = a exp — . ^ 
s ^ \ /ST 



+ h 



for T > Tc, (18) 



where b = 0.061. The parameter a controls the height and 
AT controls the width of the C/s curve and are given by 



= 0.901, AT = 



14.5' 



(19) 



We will change these values to explore the various cases 
of C/s to account for the uncertainty of the height and 
width of the curve. In FIGJ2]we show the change in bulk 
viscosity profile by varying the width of the ^/s curve by 
keeping the height intact. 



We use the lower bound of the shear viscosity to en- 
tropy density ratio known as KSS bound [20] 



ril s = l/47r 



(20) 



in our calculations. We note that the entropy density is 
obtained from the relation 



P 



T 



(21) 



In FIG IT] we plot the trace anomaly (e — 3P)/T^ and 
C/s for desired temperature range. We also plot the 
constant value of rj/s = l/47r for a comparison. It is 
clear that the non-ideal EoS deviates from the ideal case 
(e = 3P) significantly around the critical temperature. 
Around same temperature ^/s starts to dominate over 
77/s significantly. We would like to note that these re- 
sults are qualitatively in agreement with Ref . |34j . 



C. Thermal photons 

During QGP phase thermal photons are originated 
from various sources, like Compton scattering q{q)g — > 
q{q)j and annihilation processes qq — >■ (77. Recently Au- 
renche et al. showed that two loop level hremsstrahlung 
process contribution to photon production is as impor- 
tant as Compton or annihilation contributions evaluated 
up to one loop level |S2]. They also discuss a new mech- 
anism for hard photon production through the annihila- 
tion of an off-mass shell quark and an antiquark, with 



the off-mass shell quark coming from scattering with an- 
other quark or gluon. These processes in the context 
of hydrodynamics of heavy ion collisions were studied 
in Refs.[T31 E]. Until recently only the processes of 
Compton scattering and qq- annihilation were considered 
in studying the photon production rates. 

The production rate for hard {E > T) thermal pho- 
tons from equilibrated QGP evaluated to the one loop 
order using perturbartive thermal QCD based on hard 
thermal loop (HTL) resummation to account medium ef- 
fects. The Compton scattering and qq- annihilation con- 
tribution is [IJ H [S] 



E 



dN _ 1 
d^xd^p ~ 2^' 




-^/^ In 



cE 

QsT 



(22) 

where the constant c ~ 0.23 and a and Ug are the electro- 
magnetic and strong coupling constants respectively. In 
summation / denotes the flavours of the quarks and e/ 
is the electric charge of the quark in units of the charge 
of the electron. 

The rate of photon production due to Bremsstrahlung 
processes is given by[52j 



E 



dN _ 8 
d^xd^p TT 



raas 




e-^/^iJT-JL)IiE,T) 



(23) 
where Jt ~ l-H and Jl ~ 1-06 for two flavours and 
three colors of ouarksfMj. The expression for I{E,T) is 
given by 



IiE,T) 



3C(3) 



TT^E 

~6T 



E 



ln(2) + 4 Li3 (-e-l^l/T) +2^0, (-e-l^l/^) - | In (l 



-|E|/T 



(24) 



and Li are the polylogarithmic functions given by 

+00 
Lia(z) = J2 



z 
n** 



Now the rate due to qq- annihilation with an additional 
scattering in the medium is given by, 



E 



dN _ 8 
d'^xd^p 37r^ 



aas 



J2e}] ETe-""^^ [Jt 



Jl). 



We use the parametrization of cts{T) by Karschl53] 

6tt 



as{T) = 



(33 - 2Nf) ln(8T/Tc) 



(25) 
3]: 

(26) 



for our rate calculations. Here Nf is the number of 
quark flavors in consideration. 

In FigJ3) we plot the different photon rates for a flxed 
temperature T = 250 MeV. It shows the contributions 
from Bremsstrahlung (Brems), annihilation with scat- 
tering (A-l-S) and Compton scattering together with 
qq- annihilation (C-l-A). Bremsstrahlung contributes to 
the photon production rate upto E ^ 1 GeV only, 
afterwards A-(-S and C-l-A processes become dominant. 
This observation is in complete agreement with with 
Ref.[Il]. 

The total photon rate is obtained by adding different 
temperature depended photon rate expressions. Once 
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FIG. 3. Hard thermal photon rates in QGP as a function of 
energy for a fixed temperature T=250 MeV. Photon rates are 
plotted for different relevant processes. 



TABLE I. Initial conditions for RHIC 
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the evolution of temperature is known from the hydro- 
dynamical model, the total photon spectrum is obtained 
by integrating the total rate over the space time history 
of the collisionIS? , 



dN 
d?p±dy 



= / d^x 



E 



dN 
d^pd'^x 

Vm,c 



(27) 



= g / dr T 

J To 



dy [E 



Vm, 



dN 

d^pd'^x 



where tq and ti are the initial and final values of time 
we are interested, ynuc is the rapidity of the nuclei 
whereas Q is its transverse cross-section. For a Au 
nucleus Q ~ 180/to^. p± is the photon momentum in 
direction perpendicular to the collision axis. The quan- 



tity ( 



E 



dN 



is Lorentz invariant and it is evaluated 



in the local rest frame in equation (27 1. Now the photon 



energy in this frame, i.e., in the frame comoving with 
the plasma, is given as p±cosh{y — y'). So once the ra- 
pidity and p±^ are given we get the total photon spectrum. 



III. RESULTS AND DISCUSSION 

In order to understand the temporal evolution of 
temperature T{t), pressure P{t) and viscous stresses - 



$(t) and n(T), we numerically solve the hydrodynamical 
equations describing the longitudinal expansion of the 
plasma: ( 10p2 ). We use the non-ideal EoS obtained 
from equations (161 and (17). Information about vis- 



cosity coefficients C, and ry are obtained from equations 
using equation (21 1. We need to specify the 



PP ) 



initial conditions to solve the hydrodynamical equations, 
namely the initial time tq and Tq. We use the initial 
values relevant for RHIC experiment given in Table I, 
taken from Ref. [13] . We will take initial values of viscous 
contributions as $(ro) = and n(To) = 0. We would like 
to note that our hydrodynamical results are in complete 
agreement with that of Ref. |3S] . 

Once we get the temperature profile we calculate 
the photon production rates. Total photon spectrum 
-^dTWx (*^ ^ function of rapidity, y and transverse 
momentum of photon, pj_) is obtained by adding 
different photon rates using equations (22),(23),(25) 



and convoluting with the space time evolution of the 
heavy- ion collision with equation (27). The final value 



of time Ti is the time at which temperature evolves to 
critical value rj, i.e.; T(ti) — Tc. In all calculations 
we will consider the photon production in mid-rapidity 
region (y = 0) only. 

We will be exploring various values of viscosity and 
its effect on the system. Since there is an ambiguity 
regarding the height and width of C/s curve, we will 
vary the parameters a and AT from its base value given 
in equation (19). By this we will able to study the effect 
of variation of C on the system. The varied values of 
the parameters are represented by a' and AT'. We note 
that unless specified we will be using the base values 
of bulk viscosity parameters ( |19[ ) in our calculations. 
Throughout the analysis we will keep the shear viscosity 
1] to its base value given by equation ( 20 ) . 



In order to understand the effect of non-ideal EoS in 
hydrodynamical evolution and subsequent photon spec- 
tra we compare these results with that of an ideal EoS 
(e = 3P). We consider the EoS of a relativistic gas of 
massless quarks and gluons. The pressure of such a sys- 
tem is given by 



P = aT^ 



16 



21 



'-N 



90 



(28) 



where Nf = 2 in our calculations. Hydrodynamical evo- 
lution equations of such an EoS within ideal (without 
viscous effects) Bjorken flow can be solved analytically 
and the temperature dependence is given by [36] 



T^T( 



"(t) 



1/3 



(29) 



where tq and Tq are the initial time and temperature. 
While considering the viscous effect of this ideal EoS, 
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FIG. 4. Temperature profile using massless (ideal) and non- 
ideal EoS in RHIC scenario. Viscous effects are neglected 
in both cases. System evolving with non-ideal EoS takes a 
significantly larger time to reach Tc as compared to ideal EoS 



scenario. 



FIG. 5. Figure shows time evolution of temperature with non- 
ideal EoS for different combinations of bulk (11) and shear (<1>) 
viscosities. Non zero value of bulk viscosity refers to equations 



( TsflQ I and non zero shear viscosity is calculated from equa- 



tion (201 



we will solve the set of hydrodynamical equations ( |10| - 
11 1, since effect of bulk viscosity can be neglected in the 
relativistic hniit when the equation of state P — e/3 is 
obeyed [5T] . 



Hydrodynamics with non-ideal and ideal EoS 

FIG. Q shows plots of temperature versus time for the 
ideal and non-ideal equation of states. The temperature 
profiles are obtained from the hydrodynamics without 
incorporating the effect of viscosity. The figure shows 
system with non-ideal EoS takes almost the double time 
than the system with ideal massless EoS to reach Tc- 
So even when the effect of viscosity is not considered, 
inclusion of the non-ideal EoS makes significant change 
in temperature profile of the system. This can affect the 
corresponding photon production rates (below). 

Now we analyse the viscous effects. Role of shear vis- 
cosity in the boost invariant hydrodynamics of heavy ion 
collisions, for a chemically nonequilibrated system, was 
already considered in Ref.[17 . 

Next we consider possible combinations of $ and 11 in 
non-ideal EoS case and study the corresponding temper- 
ature profiles as shown in FIG. ([5|. As expected viscous 
effects is slowing down temperature evolution. For the 
case of non zero bulk and shear viscosities (11 7^ 0; $ 7^ 
0), temperature takes the longest time to reach Tc as 
indicated by the top most curve. This is ~ 1.5 fm/c 
greater than the no viscosity case (the lowest curve). The 
remaining two curves show that the bulk viscosity dom- 




FIG. 6. Longitudinal pressure Pz for various viscosity cases 
shown in FIGH 



inates over the shear viscosity when the value of T ap- 
proaches Tc and this makes the system to spend more 
time around Tc- However the intersection point of the 
two curves may vary with values of a and AT as high- 
lighted by FIGH 



Non-ideal EoS and Cavitation 

Let us note the fact that 11 < [3S]. From the defini- 
tion of longitudinal pressure P^ = P -I- 11 — $ it is clear 
that if either C (11) or 77 ($) is large enough it can drive Pz 
to negative values. Pz = defines the condition for the 
onset of cavitation. At this instant when of Pz becoming 
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FIG. 7. Temperature is plotted as a function of time. With 
peak value (a) of (/s remains same while width (AT) varies. 
Solid line in the curve ends at the time of cavitation, while 
the dashed lines shows that how system would continue till 
Tc if cavitation is ignored. Figure shows that larger the AT 
shorter the cavitation time. 



zero the expanding fluid will break apart in to fragments 
and hydrodynamic treatment looses its validity (see for 
e.g. Ref.[55]). Recent experiments at RHIC suggest r]/s 
to its smallest value ^ I/Att. And such a small value of 
77/s alone is inadequate to induce cavitation. Therefor 
we vary the bulk viscosity values by changing a and AT 
to study the cavitation. In the discussion that follows we 
will use Tc to denote the time when cavitation occurs. 

In FIGs|6]and[7]we plot P^ and T as functions of time 
for different values of AT while keeping a (=0.901) fixed. 
It shows that higher value of AT is leading to a shorter 
cavitation time. For the values of a and AT given by 
equation ( |19[ ) we find that around Tc = 2.5 fm/c, Pz 
becomes zero as shown by the curve at the bottom of 
the FIGJ6) In this case cavitation occurs when system 
temperature is larger than T^. This can be seen from 
the top curve of the in FIGJTJ End point of the solid 
line in the top curve occurs at T ^ 210 MeV and 
Tf, — 2.5 fm/c. Had we ignored the cavitation, system 
would have taken a time Tf — 5.5 fm/c to reach Tc which 
is significantly larger than Tc as seen from FIGJT] This 
shows that cavitation occurs rather abruptly without 
giving any sign in the temperature profile of the system. 
The hydrodynamic evolution without calculating P^ 
may end up in over estimating the evolution time and 
subsequent photon production. 

A similar analysis can be carried out by keeping AT 
fixed (= Tc/14.5) and varying parameter a. We show the 
cavitation times corresponding to changes in a and AT 
(denoted by a' and AT') in FIGJS) The dashed curve in 
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FIG. 8. Cavitation time t^ as a function of different values of 
height (a') and width (AT') of C/s curve. 
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FIG. 9. Photon flux as function of transverse momentum for 
different equation of states. No effect of viscosity included in 
the hydrodynamical equations. 



FIGJ8]shows Tc as a function of a, while keeping AT fixed. 
The curve shows that r^ decreases with with increasing 
a. Solid line shows how Tc varies while keeping a fixed 
and changing AT. 



Thermal Photon Production 

We have already seen that the calculation of photon 
production rates require the initial time tq, final time ri 
and T{t). ti and T{t) are determined from the hydro- 
dynamics. Generally ri is taken as the time taken by the 
system to reach Tc, i.e.; Tf. But when there is cavita- 
tion, we must set ri = Tc. Therefor photon productions 
will be influenced by cavitation, temperature profile and 
non-ideal EoS near T^. 
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FIG. 10. Photon spectrum obtained by considering the effect 
of cavitation (dashed line). For a comparison we plot the 
spectrum without incorporating the effect of cavitation (solid 
line). 



FIG. ^ shows the photon production rate calculated 
using ideal (massless) and non-ideal EoS. The figure 
shows that non-ideal EoS case can yield significantly 
larger photon flux as compared to the ideal EoS. At en- 
ergy E = 1 GeV, photon flux for the non-ideal EoS is 
60% larger than that of ideal EoS case. This is because 
the calculation of the photon flux is done by performing 
time integral over the interval between the initial time 
To and the final time Tf. tq is same for both the system 
while the Tf for the case with non-ideal EoS is two times 
larger than the ideal EoS. Since the non-ideal EoS allows 
the system to have consistently higher temperature over 
a longer period as compared to the massless ideal-gas 
EoS, more photons are produced. 

Next we try to observe the effect of cavitation in pho- 
ton production. We emphasis that rates should only be 
integrated up to r^. In FIGjTOJ photon rates are calcu- 
lated for the two cases. In the dashed curve the effect of 
cavitation is taken into account and ti = Tc = 2.5 fm/c. 
The solid line represent the same case but with the ef- 
fect of the cavitation is ignored and ti — Tf — 5.5 fm/c. 
We see from the solid curve that we end up over esti- 
mating the photon rates at p± = 0.5 GeV by ~ 200 % 
and p± — 2 GeV over estimation is about 50%. So it 
is clear that information about cavitation time is crucial 
for correctly estimating thermal photon production rate. 

In the FIG[T1] we plot photon production rates for 
various cavitation times obtained by varying AT (with 
a — 0.901 is fixed). Here the enhancement in the photon 
production when AT is reduced to half of its base value 
is ~ 75% at pA^ = 0.5 GeV and - 55% at p_L == 1 GeV. 



A further reduction of the parameter value to Ar/4 is 
enhancing the photon production by ~ 110% at p± — 0.5 
GeV and - 80% at p_L = 1 GeV. Reduction in AT 
amounts to increase in the cavitation time, which in turn 
would increase the time interval over which photon pro- 
duction is calculated. Therefor this increases the photon 
flux. 
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FIG. 11. Photon production rates showing the effect of dif- 
ferent cavitation time. 



IV. SUMMARY AND CONCLUSIONS 

Thus using the second order relativistic hydrodynamics 
we have analyzed the role of non-ideal effects near Tc 
arising due to the equation of state, bulk-viscosity and 
cavitation on the thermal photon production. Since the 
experiments at RHIC imply extremely small values for 
rj/s, the shear viscosity play a sub dominant role near T^ 
in the photon production. 

We have shown using non-ideal EoS using the recent 
lattice results that the hydrodynamical expansion gets 
significantly slow down as compared to the case with the 
massless EoS. This in turn enhances the flux of hard ther- 
mal photons. 

Bulk viscosity play a dual role in heavy-ion collisions: 
On one hand it enhances the time by which the system 
attains the critical temperature, while on the other hand 
it can make the hydrodynamical treatment invalid much 
before it reaches Tc. We have shown that if the phe- 
nomenon of cavitation is ignored one can have erroneous 
estimates of the photon production. Another result we 
would like to emphasize is that reduction in cavitation 
time can lead to significant reduction in the photon pro- 
duction. We hope that this feature may be useful in 
investigating the signature of cavitation. 
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